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[57] ABSTRACT 

A non-invasive flow analysis system and method wherein a 
sensor, such as an acoustic sensor, is coupled to a conduit for 
transmitting a signal which varies depending on the char- 
acteristics of the flow in the conduit. The signal is amplified 
and there is a filter, responsive to the sensor signal, and tuned 
to pass a narrow band of frequencies proximate the resonant 
frequency of the sensor. A demodulator generates an ampli- 
tude envelope of the filtered signal and a number of flow 
indicator quantities are calculated based on variations in 
amplitude of the amplitude envelope. A neural network, or 
its equivalent, is then used to determine the flow rate of the 
flow in the conduit based on the flow indicator quantities. 

50 Claims, 5 Drawing Sheets 
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FLOW ANALYSIS SYSTEM AND METHOD 

RELATED INVENTIONS 

This application is a continuation-in-part application of 
U.S. patent application Ser. No. 08/333,213, filed Nov. 2, 
1994, U.S. Pat. No. 5,600,073. 

GOVERNMENT RIGHTS 

This material is based upon work supported by NASA 10 
Contract Number NAS- 13-691. The Government, therefore, 
has certain rights in this invention. 

FIELD OF INVENTION 

This invention relates to a method and system for ana- 15 
lyzing different types of flows in a conduit including single 
or two phase flows and for determining the mass flow of a 
single phase flow or quality and mass flow rate of a two 
phase flow using non-invasive sensors. 

20 

BACKGROUND OF INVENTION 

A two phase flow in a conduit is a flow which includes 
both gas and liquid, or gas and solid, or liquid and solid. An 
example of a gas/liquid two phase flow is water and air 25 
flowing in a pipe; an example of a gas/solid two phase flow 
is coal particles and air flowing in a pipe. Ultrasonic 
methods for determining the presence of a two phase flow in 
a conduit are known. See, e.g., US statutory invention 
registration No. H608. Essentially, an ultrasonic pulse is sent 30 
transversely through a pipe and if the flow is single phase 
(i.e. all liquid), a return pulse is detected after a time lag as 
a return echo which is strong and reasonably sharp off the far 
wall of the pipe. If gas bubbles are entrained in the liquid, 
there are multiple small reflections and diffusion or attenu- 35 
ation of the main return echo off the far wall of the pipe. If 
a two phase flow with a defined gas/liquid interface is 
present in the pipe, the return echo is fairly strong but earlier 
in time than in the situation where there is only liquid flow 
since the return echo bounces off the gas/liquid interface ^ 
instead of the far wall of the pipe. Finally, if the flow has a 
defined gas/liquid interface and also gas bubbles in the 
liquid, there are multiple small reflections due to the bubbles 
in the liquid and the return echo is both attenuated and 
earlier in time than would be the case with only liquid flow. 45 
Such measurement methods, however, which only detect the 
presence of a two-phase flow, do not completely define the 
two phase flow. 

Quality is the mass fraction of the two-phase flow that is 
in the gaseous phase. Quality together with mass-flow 50 
determines the amount of energy (enthalpy) that is con- 
vened by the flow and thus is a key variable used to define 
the status of the flow system. Accordingly, quality and/or 
mass flow measurements are needed to fully define the flow. 
One reason that mass flow rate and the quality measure- 55 
ments are needed is to adjust the rate of one phase of the flow 
in a system. 

A typical liquid/gas two phase flow comprises a liquid 
film in contact with part or all of the conduit wall (depending 
on flow parameters and flow orientation with respect to 60 
gravity). The liquid is largely separate from a continuous or 
intermittent vapor flow. Since an identifiable liquid-vapor 
interface exists, it is possible to analyze the geometry, flow 
rate, and axial pressure drop behavior of the liquid and vapor 
flow separately, equating boundary conditions as appropd- 65 
ate. If the thickness of the liquid flow in the conduit can be 
determined, various flow models can be used to predict fill 
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thickness versus quality for a number of mass flow rates. 
See, for example. Wallis, G. B., One-dimensional Two- 
Phase Flow, McGraw-Hill, New York, N.Y., 1969. pages 
51-54, and 315-374; Lockhart, R. W. and Martinelli, R. W, 
“Proposed Correlation of Data for Isothermal Two Phase, 
Two Components Flow in Pipes”. Chemical Engineering 
Progress, Volume 45, No. 1, 1949, pages 39-48; and 
Deissler, R. G„ “Heat Transfer and Fluid Friction for Fully 
Developed Turbulent Flow of Air and Super Critical Water 
with Variable Fluid Properties”, Transactions, ASME vol- 
ume 76, No. 1, 1954. page 73. 

But the precursor step of detecting the thickness of the 
liquid flow using ultrasonic methods is troublesome. The 
presence of bubbles of gas in the liquid flow, the presence of 
large waves of liquid traveling in the conduit, small scale 
thickness changes in the liquid/vapor interface, and other 
similar “chaotic” conditions within the conduit severely 
affect the ability to determine fill thickness using ultrasonic 
techniques. If one or more of these conditions are present 
within the conduit, a plot of the return echoes from an 
ultrasonic transducer is not a good indicator of film thick- 
ness. Moreover, a low flow rate with a high quality results 
in a highly chaotic flow as does a high flow rate with a low 
quality. Such chaotic flows render known film thickness 
measurement techniques unreliable. 

Therefore, a trace of the return echoes from such a chaotic 
flow alone is seemingly not a good indicator of film thick- 
ness. Other techniques for measuring the thickness of the 
liquid film that is usually in contact with the wall of the 
conduit include sampling, thermal probes, film conductivity 
or capacitance measurements, and gamma densitometry. 
Although each of these techniques exhibits strengths and 
weaknesses, no technique offers the advantages of 
reflective-mode ultrasound. Ultrasound techniques are non- 
invasive, offer rapid response, excellent long term accuracy 
and sensitivity, and are applicable to all working fluids over 
a very broad range of temperatures. Moreover, even if 
non-ultrasonic thickness measurement techniques are used, 
the various flow models used to evaluate flow quality and 
mass flow rates are based on a number of assumptions which 
can lead to inaccuracies. On the other hand, quality and/or 
mass flow measurements cannot be accurately taken without 
measuring film thickness or a related parameter, void frac- 
tion. Flow meters, for example, do not indicate how much of 
the flow is liquid or gas and flow meters cannot be used in 
all situations. 

A typical gas/solid two phase flow, such as coal particles 
entrained in an air flow, generally comprises a rope like 
structure of coal particles travelling in the pipe. There are no 
current techniques which accurately measure the mount of 
coal in the pipe. Trial and error methods commonly used in 
coal power plant operation, can result in poor efficiency and 
air pollution. In order to optimize combustion, the amount of 
coal and the amount of air delivered to the burner must be 
known. 

Therefore, in addition to determining the presence of a 
two phase flow, considerable research has been performed 
on various means of actually measuring two-phase flows. 
These efforts have largely attempted to characterize an 
average value of some aspect of the flow, such as a pressure 
drop, void fraction, fill thickness, velocity, or density. One 
problem with this approach is that knowledge of any single 
value is not sufficient to define a two-phase flow. TXvo-phase 
flows comprise two separate flows (of phase A and Phase B) 
that interact in extremely complex ways. If average values 
are used, at least two independent quantities must be mea- 
sured to define the flow. In addition, a given pair of 



5 , 741,980 


3 

observations, such as a pressure drop and a velocity, often 
does not provide a sensitive indication of the flow rates of 
phase A and phase B for a broad range of conditions. Thus, 
different combinations of observations are often needed for 
different flow conditions. 

There are applications that can be well served by suitably 
developed instruments based on currently known averaging 
techniques. However, there are many more applications for 
which combinations of currently available averaging mea- 
surements will not provide desirable results. For example, 
some applications demand a completely non-in vasive flow 
measurement Others may be geometrically constrained, so 
that only instruments of a given size or shape may be used. 
Other applications may require accurate flow measurement 
over an extremely broad range of flow conditions. Still 
others may be very cost-driven, so that the instrument must 
be very inexpensive. The current invention offers the advan- 
tage that any meaningful measurement technique, used 
rapidly and persistently, can be used to determine the flows 
of both phases. Since every application permits at least some 
meaningful flow observation to be made, the current inven- 
tion ensures that a practical instrument can be developed. 

The approach of the subject invention arises from simple, 
but profound observations about two-phase flows. First, they 
are deterministic, in that they satisfy the laws of physics. 
Thus, while their evolution is extremely complex, there is an 
underlying order to the flow behaviors. The behavior of a 
given wave, particle, or bubble, although complicated, is not 
truly random. 

Second, since they are deterministic and behave in a 
complex fashion, they are likely to be chaotic. The word 
“likely” is used here because it has not yet been proven in 
generality that fluid turbulence satisfies the mathematical 
definition of chaos. It is not yet known what (if any) kinds 
of complex behaviors are possible that are neither chaotic 
nor random. It currently is not possible to make generalities 
about the behaviors of such a system. For want of a more 
conclusive answer from topologists, and since the view of 
two-phase flows as chaotic is consistent with the evidence so 
far available, the current invention assumed that the two- 
phase flows of interest are chaotic. 

Two-phase flows are dissipative, i.e. given the 
opportunity, any work that is imposed on them is eventually 
lost to viscosity. Thus, given an arbitrary initial condition 
(e.g., means of mixing at the inlet to the flow conduit), a 
two-phase flow will settle into a pattern of behavior that is 
similar to that of other flows with the same flows of phase 
A and phase B but different initial conditions. Actually, it is 
difficult to prove this conclusively without generating many 
flow conditions with a variety of initial conditions and 
comparing their properties in detail. However, it is a basic 
tenet of the arts of fluid mechanics and two-phase flow that 
this is the case: if this were not the case, it would not be 
possible to generate models or correlations of two-phase 
flow behaviors. This is also consistent with the properties of 
dissipative chaotic systems. 

If two-phase flows are dissipative and chaotic, then key 
statements can be made about their behaviors. Principal 
among these is the existence of a single underlying behavior, 
a “strange attractor”, the shape of which changes as flow 
parameters change. The strange attractor is an extremely 
complicated path (in mathematical phase space) that defines 
all trajectories of the system with time. It is limited to a finite 
portion of phase space and is a single unending, open (i.e., 
never repeating) path in which points that are initially near 
one another diverge rapidly from one another with time 


4 

(called sensitivity to initial conditions). The conclusions that 
can be drawn from the existence of a strange attractor are 
far-reaching. The principal conclusion for current purposes 
is that any observation of the system behavior, made over a 
5 period of time, is a mapping of the strange attractor. If the 
observation is made with constant time increment between 
measurements, it is a smooth mapping. Any smooth map- 
ping of a strange attractor contains an mount of information 
about the system behavior that is comparable to any other 
10 smooth mapping with the same time increment and mea- 
surement sensitivity. Thus, any of a variety of measurement 
methods may be used with a two-phase flow with equal 
conviction that meaningful information is obtained. 

This line of discussion is fairly well established in the art 
15 of chaos theory, but by itself is not sufficient to permit the 
measurement of two-phase flows. The reason for this is that 
the argument does not disclose how the evolution of the flow 
observations can be related to the flow conditions. In fact, 
various researchers have attempted to relate two-phase flow 
20 conditions to time-series measurement. The best known of 
these are Jones, O. C. and Zuber, N., ‘The Interrelation 
between void fraction fluctuations and flow patterns in 
two-phase flow” Int’l J. of Multiphase Flow, v2. page 
273-306, 1975; Hubbard, M G., and Dukler, A. E., ‘The 
25 Characterization of Flow regimes for Horizontal two-phase 
flow: 1 . statistical analysis of wall pressure fluctuations”. 
Proceedings of the 1966 Heat Transfer and Fluid Mechanics 
Institute, Saad, M. A. and Miller, J. A. eds.. Standard 
University Press, pages 100-121, 1966. Jones and Zuber 
30 identified liquid-gas flow regimes from the probability den- 
sity function of X-ray attenuation measurements. Hubbard 
and Dukler identified flow regimes from frequency spectra 
of pressure signals from liquid-gas two-phase flows. In 
neither of these cases was the flow rate of either or both 
35 phases determined. 

The technical literature has many references describing 
efforts to examine or develop instruments of various kinds to 
measure or characterize two-phase flows. Overviews include 
Hsu, Y. Y., and Graham, R. W., Chapter 12: Instrumentation 
40 for Two-Phase Flow, in Transport Processes in Boiling and 
Two- Phase Systems ; McGraw-Hill, 1976; and Mayinger, F., 
Chapter 16: Advanced Optical Instrumentation, in Two- 
Phase Flow and Heat Transfer in the Power and Process 
Industries , Bergles. A. E., et al, editors. Hemisphere Pub- 
45 lishing Company, 1981. The bulk of these efforts have 
generated results of sufficiently limited application to have 
remained largely of research interest True two-phase 
flowmeters, i.e., instruments that purport to define the flows 
of both phases A and B, are not widely available on the 
50 commercial market. 

In industrial practice, the most widely used instruments 
for two-phase flows have been photon attenuation instru- 
ments. These instruments determine the attenuation of pho- 
tons (typically microwaves or gamma rays) as they pass 
55 through the flow. The greater attenuation of one phase than 
the other (essentially from higher density) is used to char- 
acterize the portion of the flow channel cross section that is 
filled with each phase. Alternatively, this may be viewed as 
characterizing the average density of the flow. Depending 
60 upon the specific geometry of the instrument, it may be 
rendered more or less sensitive or the distribution of the 
phases across the flow channel, and thus may be used to 
identify the flow regime (e.g., bubbly, slug, stratified, 
annular, or mist flows). These instruments provide only a 
65 rough indication of the flow condition, because their sensi- 
tivity to the amount of each phase that is present is limited. 
In particular, gamma densitometers are highly sensitive to 
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even trace quantities of lead in a flow, severely limiting their 
accuracy in many applications of interest (most notably in 
petroleum pipelines). Despite the limitations of these instru- 
ments in their originally intended embodiment, they could 
be used to advantage with the current invention to accurately 5 
determine the flows of phases A and B. 

To generate useful information about the flow rates of 
phases A and B, an attenuation measurement must be 
combined with some indication of flow velocity. Such an 
indication may be obtained by making attenuation measure- i0 
ments at two closely spaced stations along the flow duct and 
cross-correlating the resulting signals. The time delay of the 
peak in the cross-correlation curve corresponds to an 
approximate time delay for flow propagation. Dividing the 
spacing of the instruments by this time delay provides a 15 
characteristic velocity. The average flow density and this 
velocity can be correlated to the flow rates of phases A and 
B. Correlation is needed to correct for the inevitable “slip” 
that occurs between phases (because they do not flow with 
identical velocity). 20 

Even with suitable calibration, the accuracy of cross- 
correlated attenuation measurements is limited because of 
the poor sensitivity of the density measurement for many 
flow conditions. This limited accuracy is implied in U.S. Pat. 

No. 4,683.759 to Skarvaag et al, wherein, this basic idea is 2 5 
used to measure liquid-gas two-phase flows. However, the 
determination of liquid and gas flow rates is discussed for 
only one specific flow regime, called slug flow, in which the 
liquid and gas flows are largely intermittent, and the peak in 
the cross-correlation function quite sharp. 30 

Other instrument systems have been devised to observe 
two-phase flows. In U.S. Invention registration H608 to 
Goolsby, an ultrasonic measurement technique is used to 
determine whether gas is present in a liquid flow. In this 
instrument, echo-mode ultrasound is used to determine the 35 
location of the second major reflection interface (the first 
being between the liquid and the tube wall). If the time-off- 
flight of the acoustic wave is less than that associated with 
a full tube (second reflection from the far wall), then a 
liquid-gas interface is present. Actually, this approach has 40 
been used to study two-phase flows quantitatively for quite 
a few years. 

In U.S. Pat. No. 4,193,291 to Lynnworth, an ultrasonic 
method of determining flow density is described. This 
technique is based on the different attenuation rates of 45 
torsion waves in a body depending upon the density of the 
fluid in which the body is immersed. Various embodiments 
are described that render the instrument more or less sensi- 
tive to the distribution of the phases in the flow duct. This 
instrument is limited to liquid — liquid or liquid- gas two- 50 
phase flows. One unfortunate aspect of this instrument is its 
intrusiveness into the flow. The protrusion of the instrument 
into the flow raises the potential for damage to the instru- 
ment from debris carried by the flow, generates an undesir- 
able pressure drop and flow disruption, and requires seals 55 
between the instrument and the flow duct which reduce the 
reliability of the flow system. Another unfortunate aspect of 
this instrument is that no single embodiment is described 
that determines the average density across the entire flow 
cross section. Thus, each embodiment is applicable to a 60 
limited range of flow regimes. Further, the issue of fluid 
wetting is not addressed: if the liquid wets the material of the 
sensor, the apparent density may be skewed strongly toward 
the liquid density. Even with these limitations, the mecha- 
nism of this measurement approach could be used with the 65 
current invention to provide an accurate determination of the 
flow rates of phases A and B. 
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The aforementioned patent states that the density mea- 
surement may be combined with acoustic velocimetry to 
determine the flow rates of both phases. In acoustic 
velocimetry, an acoustic wave is propagated by one trans- 
ducer downstream through the flow and its time-of-flight to 
another transducer measured. A second wave is propagated 
upstream to determine the propagation time against the flow. 
Comparison of these propagation times determines both the 
effective speed of sound of the flow and its propagation 
velocity. If acoustic velocimetry works in twophase flows at 
all, the measurements would be very sensitive to the flow 
regime. For example, in annular flow, a continuous liquid 
film is in contact with the flow duct wall. Any acoustic 
waves that enter the flow will effectively “short circuit” 
through the liquid (with its very high sonic velocity and 
relatively low attenuation) so that only the liquid velocity 
(skewed by the acoustic wave path) would be measured. By 
contrast, in stratified flow, a liquid flow on the bottom of the 
flow duct is effectively separated from a gas flow in the top 
of the duct. If acoustic waves traveling through the liquid 
can couple sufficiently with the gas flow, then a velocity that 
is an average of the liquid and gas velocities will be 
measured. This velocity would be quite different from (much 
higher than) that of an annular flow, even though annular 
flows often occur at much higher velocities. Thus, the 
calibration of the instrument output with different flow 
conditions would involve some very strong nonlinearities 
likely to result in poor accuracy. 

U.S. Pat. No. 4,991,124 to Kline describes a different 
ultrasonic instrument for determining a fluid density. This 
technique is based on determining the velocity of sound and 
rate of attenuation of acoustic energy in the fluid. Because 
this technique relies on multiple reflections of the acoustic 
energy, which would be extremely difficult to detect in 
two-phase flows, it probably could not be applied to a 
two-phase flow. 

AEA Technology, of the United Kingdom, has publicized 
a two-phase flowmeter for use in oil and gas fields 
[Anonymous, “Non -Intrusive Meter Measures Oil and Gas 
Flows”, Competitive Edge, Issue 4, pg. 3, Spring, 19941. 
This instrument uses a pulsed neutron beam which counts 
hydrogen, carbon, oxygen, and chlorine atoms passing the 
sensing point Short bursts of radiation are used to activate 
oxygen atoms, which can be tracked as they move to define 
a flow velocity (a second measurement method). This instru- 
ment employs two averaging techniques to determine the 
flow rates of (potentially) several phases. However, it 
depends upon the phases being of distinct compositions to 
define their flow rates separately. If the two phases were of 
the same composition (a so-called single-component two- 
phase flow), then only a total flow measurement would be 
obtained. While this system may prove effective for its 
intended application in oil fields, its cost, complexity, and 
operational limitations will limit its use elsewhere. 

Finally, the analysis of liquid and two-phase flows of 
cryogenic propellants presents a number of unique prob- 
lems. The fact that the flow is unavoidably two-phase during 
start up and in certain transients is one major problem 
Intrusive mechanical flow meters can not accurately analyze 
two-phase flows. Turbine meter moving parts can be 
destroyed, and even orifice plates and venturies can suffer 
significant erosion. These problems are particularly severe 
in the case of liquid oxygen flows which can rapidly destroy 
welded components. Intrusive flow measurement devices 
also require the minimization of the pressure drop and can. 
in some situations, cause a severe hazard. For example, 
coriolis mass flow meters are notoriously unreliable with 
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two-phase flows. See Baird, R. S., “Flow Meter Evaluation 
For On-Orbit Operations,” NASA technical memorandum 
100465, August 1988. Magnetic flow meters do not work 
with conductive fluids or two-phase flows. Also, past efforts 
to use triboelectric measurements have also proven to be 5 
unsatisfactory. See Bematowicz, H.; Cunningham, J. and 
Wolff. S., “Mass Flow Meter Using The Triboelectric Effect 
For Measurement In Cryogenics,” NASA CR- 179572, April 
1987; and Dechene, R„ “Mass Flow Measurement Of Liq- 
uid Cryogens Using The Triboelectric Effect.” NASA 
CR- 1795 19. August 1986. One current effort to use active 
ultrasonic measurements suffers the limitation of signal loss 
during the inverted angular flow transient that is common 
with system start-up, and the loss of accuracy with the 
significant gaseous void fraction. 

Accordingly, measurement of liquid and two-phase flows 
of cryogenic propellants and other fluids in the oil and gas 
industry, the chemical processing industry, and the nuclear 
and fossil fuel power generation industry requires a non- 
intrusive sensor, operable at temperatures ranging from the 
cryogenic level to high levels which provides a reliable 
measurement of both liquid and two-phase flows. 

Unfortunately, no such system exists at the present time. 

SUMMARY OF INVENTION 

It is therefore an object of this invention to provide a 
highly reliable and efficient flow analysis system and 
method. 

It is a further object of this invention to provide such a 
flow analysis system and method which is non-intrusive. 

It is a further object of this invention to provide such a 
flow analysis system and method which is capable of 
operating at cryogenic or elevated temperatures and in 
adverse conditions. 

It is a further object of this invention to provide such a 
flow analysis system and method which provides reliable 
measurement of both liquid and two-phase flows. 

This invention results from the realization that the flow 
rate, the flow quality, and other flow parameters of a single 
phase or a two phase flow in a conduit can be accurately and 
efficiently determined by using a passive sensor such as an 
acoustic sensor or a strain gauge type transducer, filtering the 
sensor output to pass only the resonant behavior of the 
sensor, capturing the amplitude envelope of the filtered 
signal, sampling the amplitude envelope, calculating a num- 
ber of flow indicator quantities based on variations in the 
amplitude of the signal, and using a trained neural network 
to characterize the flow based on the flow indicator quanti- 
ties. 

This invention features a flow analysis system and 
method. The system includes a sensor in communication 
with a conduit which has a flow therein. The sensor transmits 
a signal which varies depending on the characteristics of the 
flow in the conduit. There is a filter, responsive to the sensor 
signal, and tuned to pass a narrow band of frequencies 
proximate to the resonant frequency of the sensor. There are 
means for generating an amplitude envelope representative 
of the sensor signal, and there are also means, responsive to 
the amplitude envelope, for calculating a plurality of flow 
indicator quantities representative of the flow in the conduit 
Such means may include the computer program included 
herein. Finally, there are means, responsive to the flow 
indicator quantities, for determining the flow rate of the flow 
in the conduit . Preferably, such means include a trained 
neural network which receives the flow indicator quantities 
as input and outputs, for example, the liquid flow rate and 
the gas flow rate of a two-phase flow in a conduit. 
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The means for calculating preferably includes means for 
sampling the amplitude envelope at a rate less than or equal 
to the resonant frequency of the sensor to efficiently analyze 
the sensor signal . Typically, variations in the amplitude of 
the amplitude envelope are used to calculate the flow 
indicator quantities. One such flow indicator quantity 
includes the average value of the amplitude of the amplitude 
envelope. Another flow indicator quantity includes the stan- 
dard deviation of the amplitude. Still another flow indicator 
quantity includes the characteristic auto-correlation time of 
the amplitude envelope. Numerous other flow indicator 
quantities are calculated based on variations in the amplitude 
of the sensor signal. In a preferred embodiment, the sensor 
is an acoustic sensor. 

This invention also features a flow analysis method. A 
sensor is coupled to a conduit in order to transmit a signal 
which varies depending on the characteristics of the flow in 
the conduit A filter, responsive to the signal and tuned to 
pass a narrow band of frequencies proximate the resonant 
frequency of the sensor, is used. This filtered signal is 
demodulated and amplified and used to calculate a plurality 
of flow indicator quantities representative of the flow in the 
conduit. The flow rate of the flow in the conduit is then 
determined using the flow indicator quantities. 

DISCLOSURE OF PREFERRED EMBODIMENT 

Other objects, features and advantages will occur to those 
skilled in the art from the following description of a pre- 
ferred embodiment and the accompanying drawings, in 
which: 

FIG. 1 is a schematic view of the flow analysis of this 
invention; 

FIG. 2 is a graph depicting the variation in the amplitude 
of the sensor signal generated by the flow analysis system 
shown in FIG. 1; 

FIG. 3 is a graph showing a sampling rate of 100 kHz of 
the signal shown in FIG. 2 ; 

FIG. 4 is a graph showing the demodulated signal of FIG. 
2 and a much lower and more efficient sampling rate; 

FIG. 5 is a block diagram of the primary components of 
the flow analysis system of this invention; 

FIG. 6 is a graph depicting one particular flow indicator 
quantity used to characterize a flow in accordance with the 
subject invention; 

FIG. 7 is a graph depicting another flow indicator quantity 
used to characterize a flow in accordance with the subject 
invention; and 

FIG. 8 is a schematic view of a neural network used to 
determine the flow rate of a flow in accordance with the 
subject invention. 

Flow analysis system 10. FIG. 1, includes sensor 12 
attached to conduit 14 which includes a flow therein in the 
direction of vector 16. In FIG. 1, the flow is shown as a 
two-phase flow with liquid 18 and gas 20 components 
although this is not a necessary limitation of the subject 
invention. Sensor 12 is preferably a passive acoustic “Clam- 
pOn 2000 Particle Monitor” available from ClampOn AS of 
Norway. As an alternative, an Endevco Model 725 9 A accel- 
erometer may be used. The voltage signal from sensor 12, 
which varies as a function of the condition of the flow in 
conduit 14. is sent to amplifier 22 and then passed to narrow 
band pass filter 24 which is preferably tuned to pass a narrow 
band of frequencies proximate the resonant frequency of the 
sensor (e.g., approximately 83 kHz for the “ClampOn” 
sensor; approximately 93 kHz for the Endevco sensor). Also 
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included within flow analysis system 10 are some means for 
generating an amplitude envelope of the filtered signal, for 
example, demodulator 26. A demodulator circuit such as a 
Signetics ULN2209 and MC1496K integrated circuits may 
be used. In response to the demodulated signal, computer 5 
subsystem 28 digitizes the demodulated, filtered, and ampli- 
fied signal and calculates a plurality of flow indicator 
quantities representative of the flow in the conduit. Com- 
puter sub-system 28 includes a 16 bit data acquisition board 
capable of 100 kHz sampling. There are also some means, 
responsive to the flow indicator quantities, for determining 
the flow rate of the flow in the conduit. In a preferred 
embodiment, such means include a trained neural network 
discussed below with reference to FIG. 8. 

“ClampOn” sensor 12 has a resonant frequency of 15 
approximately 83 kHz. Its passive nature, with no power 
dissipated at the sensor location, and its piezoelectric oper- 
ating principle, offers both ruggedness and capability of 
operation in a very broad range of temperatures and is thus 
ideal for cryogenic applications. Sensor 12 produces a very 20 
low level voltage output in response to turbulence-induced 
acoustic events that occur in the flow as it passes by. Because 
filter 24 discards all behaviors that do not excite the transient 
resonant response of sensor 12, the output of filter 24 is 
fairly insensitive to structure-born mechanical vibration and 25 
audible noise. Other sensors, however, such as the ultrasonic 
sensor shown in the parent application incorporated herein 
by this reference, may also be used. 

High gain amplifier 22 boosts the low level output of the 
sensor to a voltage level signal useful for data acquisition 30 
purposes. Narrow band pass filter 24 is centered around the 
natural frequency of sensor 12. Thus, the voltage signal 
provided to demodulator 26 also has a frequency of 83 kHz 
as shown in FIG. 2. The reason for narrow band pass filter 
24 is as follows. Any linear mechanical element (such as 35 
acoustic sensor 12) has two responses to a periodic driving 
force: a steady-state solution at the frequency of the driving 
force and a transient solution at the resonant frequency of the 
mechanical element. The response to any combination of 
driving forces can be viewed as a sum of these behaviors. In 40 
a flow situation, the driving forces of interest are inherently 
transient, because they involve acoustic events resulting 
from passage of turbulence induced flow fluctuations. Thus, 
the principal flow information is almost certainly contained 
in the part of this signal that occurs at the resonant fre- 45 
quency. Spurious effects, whether from electromagnetic 
interference, amplifier noise, external mechanical noise, and 
the like occur over a broad range of frequencies. These 
behaviors, which are likely to be confusing to the analysis, 
are largely ignored by the use of band pass filter 24. 50 

The reason for demodulator 26 is explained as follows. 
The sensor signal at the output of band pass filter 24 has a 
varying amplitude and a very narrow frequency range cen- 
tered about 83 kHz. Thus, a sampling rate of at least 252 kHz 
representing four samples per period would be needed to 55 
characterize the variation in amplitude accurately. Sampling 
at a more reasonable rate, for example 100 kHz, would only 
obtain approximately 1.2 samples per period as shown in 
FIG. 3 wherein the black dots represent samples. As shown 
in FIG. 3, any sampling rate of less than 252 kHz would not 60 
capture all of the signals character. Note that the amplitude 
rise as shown at 30 would be missed as would the amplitude 
rise shown at 32. 

The subject invention includes means for calculating a 
variety of statistics that characterize different behaviors of 65 
the time variation of the signal shown in FIGS. 2 and 3. A 
few of these ‘How indicator quantities” are fairly straight 
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forward. Examples include the signal’s average value or its 
standard deviation. Other flow indicator quantities, however, 
are more involved and are intended to quantify the specific 
aspects of the signal such as the characteristic sizes, 
durations, or sequences of behaviors in the signal dynamics. 
Any sampling rate which misses the amplitude rises as 
shown at 30 and 32 could result in erroneous readings. 

A sampling rate of 252 kHz would correctly sample the 
amplitude rises at 30 and 32. Such a high sampling rate 
however, would require an extremely sophisticated data 
acquisition processing system. Same is true for a sampling 
rate of even 100 kHz. 

Accordingly, in the subject invention, demodulator 26 is 
used to generate signal envelope 40, FIG. 4. which contains 
all of the key information about the sinusoidal variations of 
the sensor signal. Accordingly, a much lower sampling rate 
can then be used to fully characterize the signal variation. 
Thus, in the subject invention, even a 30 kHz sampling rate 
can be used to capture the complete character of the signal. 
This lower sampling rate is far easier to achieve in real time 
than the higher rates required to fully capture the detail 
dynamics of a 83 kHz signal. 

In this invention, then, a liquid flow rate of a two-phase 
flow is displayed via monitor 50, FIG. 5, at 52 and the gas 
flow rate as shown at 54 may be provided in real time. 
“ClampOn” sensor 12 is clamped to conduit 14 and the 
resulting signal is sent to signal conditioning, data 
acquisition, and buffering electronics 56 which include 
amplifier 22, narrow band pass filter 34. and demodulator 
26, FIG. 1. Computer software, described in more detail 
below, operating on computer sub- system 28, FIG. 1. 
receives the resulting signal, digitizes it, and calculates a 
plurality of flow indicator quantities schematically shown in 
block 58, FIG. 5. These flow indicator quantities are then 
presented to the neural network correlation sub- system 60. 
also operating on computer sub- system 28, FIG. 1 to provide 
real time output to monitor 50, FIG. 5. 

Two seconds or more of data for each test condition may 
be collected resulting in 200,000 readings in each data file. 
Such a data collection is best for measuring flows of pure 
liquid, where the dynamics are driven by small scale turbu- 
lent eddies. Two-phase flow conditions often require a 
longer observation period because they are characterized by 
the passage of longer events such as slugs or waves of fluid. 
Stared in ASCH format, each file would consume approxi- 
mately 1.5 MB of disk space. In 241 tests, over 350 MB of 
raw data is collected. 

One such flow indicator quantity, (Avg) is the average 
value of the amplitude of the signal shown in FIG. 4. A 
second flow indicator quantity is the standard deviation of 
the signal (StDev) which reflects the harmonic power of the 
signal and is an important variable for any irregularly 
varying signal that carries significant information. A third 
flow indicator quantity is the characteristic auto-correlation 
time (Auto). This is the time scale during which the signal 
loses linear correlation and is a measure of the rate of mixing 
in a process. The auto-correlation quantity is the integral 
over all time of the signal multiplied times a timed delayed 
version of the signal. For long time delays, the auto- 
correlation approaches 

nF (1) 

where n is the number of readings in the record. For no time 
delay, the auto-correlation equals: 

lx 2 


( 2 ) 



5,741.980 


11 


12 


measured in sampling intervals, the auto-correlation quan- 
tity is: 




Hx 2 — rur 
Ij? - JLxiKM 



A fourth flow indicator quantity is the average difference, 
more specifically the average absolute value difference in 
sequential readings. This flow indicator quantity provides a 
first order indication of how rapidly the data change on the 
smallest time scale. io 

A fifth flow indicator quantity is the root mean square 
difference (RMSDif) in sequential readings, a measure simi- 
lar to the average difference, but tends to weigh larger 
changes in the sequential readings more strongly than 
smaller changes. 15 

Another flow indicator quantity is the period of average 
crossing waves variable, (APd), which is reflective of the 
average number of readings that occur between upward 
crossings of the average value of the amplitude as shown at 
60, FIG. 6. and characterized as the time length of first order 10 
signal variations. 

Another flow indicator quantity, (ARMSPd), is the root 
mean square value of the period of average crossing waves 
variable which emphasizes longer-period waves. A related 
flow indicator quantity, (ADur), is the duration of average 


crossing waves which is a measure of the average number of 
readings that occur between an upward and downward 
crossing of the average amplitude value as shown at 60, FIG. 
6. A related flow indicator quantity is the root mean square 
value of the duration of average crossing waves, 
(ARMSDur). 

The period of standard deviation crossing waves flow 
indicator quantity, (sPd), represented in FIG. 7, is a measure 
similar to the average period of average crossing waves, but 
in this case the bar is raised somewhat to look at the period 
of waves that cross the average plus the standard deviation 
as shown at 70, FIG. 7. A related flow indicator quantity is 
the root mean square value of the period of standard devia- 
tion crossing waves value, (sRMSPd). Another flow indica- 
tor quantity, (sDur), is the duration of standard deviation 
crossing waves as shown in 72. FIG. 7. Finally, a related 
flow indicator quantity is the root mean square duration of 
standard deviation crossing waves. (sRMSDur). Other flow 
indicator quantities, however, are within the scope of this 
invention. 

One version of a Quick Basic 4.5 programming language 
program running under the DOS operating system, operating 
on computer sub- system 28, FIG. 1, for calculating the 
above flow indicator quantities is included herein as follows: 
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' Program FLowiNDI: accepts coonand line input of a DataFile$ ***** or the 
‘ name of a tile containing a xiet of data tile* (PileUstS) . It opens 
1 each D&t&FileS, and calculates flow indicator qoaatitiee . The results 

• are placed in an output file naiad using the input file name plus the 

» extension ■ -OUT" . For * FileUst$ input file, the output fil* contains 
1 one line for each input file &*ned in the list. 

» Usage: 

• FL0HXUDX.m lDataFUs$l l-L meList$] 

• Without arguments, program provides a description of usage syntax. 

■ with DataFtle$ provided, program analytes the single till and puts results 

• in a new file named the same as DetaFile$ but with the extension .001. 

■ with -L FileListS provided, program analyzes each file listed in PileListS 

' and puts results in a new file with the same nams as Filddat$ but with 

• the extension .OUT. 


* Kata command that uses m t inr un available conventional wemory. 

* $ DYNAMIC 


1 Default variable type is integer in this nodul e. 

DErTNT A-Z 

1 Declare the Cotnline subprogram (which ratrisves the pro g ram command line), 

» as well as the number and type of its parameters. ***4? subprogram is 
1 provided with QuickBasic 4. 5 as an example progras- 
DECLARE SUB Comline CH. a$(). Hex} 

* Dimension arrays 

DIM aS(l TO 15) ' Contains the coouand line parameters 
DIM f${l TO 100) 1 Contains lUt of input files 

DIM array ( -30000 TO 30000) 1 Contains array of integer readings, good for 

1 2 se co nd s of data at 30kHz. 

* Gee what was typed on the command line. 

CALL Comliue IN, t$0» «) 

* Figure out whether co m n and line is workable: 

IF N * 0 GOTO Bomb ' No coasaand line arguments 
IF N > 2 GOTO Bomb r Too many arguments 

IF H • 2 AND ta$ (l) <> *-L* AID a$U) <> *.!*) goto Bomb * Wrong syntax for 
1 FilehistS operation 

GOTO NoBomb 1 No problems so far.*. 

Bomb: * Give user program usage syntax pointers, then quit. 

PRINT * 

PRINT " FLOWXNDIi Calculates signature quantities for 1 or wore 

PRINT - input data files. 

PRINT • 

PRINT • osages FLOWINDI tf>ataFil»1 l‘b ListFile) 

PRINT * 

PRINT * 1) DacaFile - path and name of single input file, must be a 

PRINT * single Column of normal integers (-32766 to +32767) . 

PRINT * 2) ListFile - path and name of text file containing a list of 

PRINT ■ paths ♦ names of files to be analyzed (single column) . 

PRINT " 3) Command lina can contain DataFlle or ListFile (not both) . 
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PRINT " 4) Results are placed in an output file named using input file * 

PRINT ■ suns with extension replaced by .OUT. * 

PRINT • • 

END 


NoBomb: 1 OK, so far *0 good. Next , open the Input list file (if there 
* is one) and create an output file. 

IF N - 1 THEN ' Only one input file. 
f$U> - a$ (1) 
iPilss * 1 


ELSE ' M- 2 . FileUetS input, gtt 11 «c o£ tileo 

FileLi*e$ - »$<2) 

ON ERROR GOTO Haadlel 
flag « 1 

OPEN FileLiat$ FOR INPUT AS Ml 
ON ERROR GOTO 0 

if lies - 0 1 Number of file* in Uat 

DO WHILE NOT BOFU) ’ Read list of input files 
INPUT #1, ffSUFilee + l) 
lFiles - IF ilea * X 

LOOP 

CLOSE #1 * Dona with PileLiat$ 

flag ■ 2 

IF iFiles « 0 THEN GOTO Handle 1 * No file names 
ENO IF 

•Next, create output file with name baaed on FlleLisc$ or DacaFile$, as 
• appropriate. 

IF » - 1 THEN 

file$ - a$ (1) 

ELSE 

filet - FileLiat$ 

END IF 

' Find extension in file name 
FOR ichar - 1 TO LEN{file$) 

IF MID$ (fil*$, ichar, l) . *.» THEN 

* everything to the right of this character roust be extension 
oFileS * LEFT$(flile$, ichar) ♦ *0171* 

END IF 
NEXT ichar 

• Plug hole for case where file didn't have an extension initially 
IF RIGKT$ (oFile$, 4) <> • .OUT" THEN oFile$ • oFile$ * •.OUT" 

1 Open output file 
ON ERROR GOTO Handle! 
flag - 3 

OPEN OFileS FOR OUTPUT AS #2 
ON ERROR GOTO 0 

* Write header line in output file 

PRINT M2, " FlleName Average StDev Auto AvgOif RMSDif 
PRINT M2, ■ APd ARHSPd * ADur ARMSDux SPd SRMSPd 
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PRINT *2, ■ S&ar SBMSDur ■ 

• Next, analyze each input file in turn and print results in oFile$. 

FOR flleNum - 1 TO i Files ♦ for e ac h date tile 

QM ERROR goto Handle2 • tile span error 
OPEN fStfilaBum) FOR. INPUT AS #1^ 

ON ERROR GOTO 0 

PRINT "Analyzing m i f${fUeNuo); ■.« 

' First, go through data act to Input data array and calculate preliminary 
1 Clow Indicators. 

Jt * 0 ' Runber o£ values found in tile (so tar) 

Sums* •* o * Accumulator of sue ot raw values 

Sum2# « 0 * Acc umul ator of sue of values squared 

SumDif# * 0 1 Accumulator for sun of absolute value changes in aeq rdga. 

SumDif2# - 0 * Accumulator for sum of squares of differences in seq rdga. 

Corr# - 0 * Accumulator for product ot current value and preceding value. 

DO WHILE NOT SOF(l) 

INPUT #1, array (JR 1 - 30000) 

Jfc - Jfc ♦ 1 

IF Jt > 1 THEN 'Not first data point 

SumDifS - SincDif# 4* ABS (value - oldvaluefc) 

SumDifl# - SunCif2# ♦ (value - olftvalusfc) * 2 
Corr# ■ Corr# + value * oldvaluefc 
END IF 

Sum St - Surat ♦ value 
Sura2# - Subi 2# > value ~ 2 
oldvaluat ■ value 

LOOP 

CLOSE #1 ' Done with raw data, so close file, 

PRINT -There are JR,* 1 points in the data set.* 

'Calculate a tew results 

Ave! « Sums / JR * Average value 

SDev! « SQR{Sum2# / J* - Ave! * 2) * Standard deviation 

Auto* - (Sum2# - J* • Avei * 2) / <Suo2# - Corx#> * Characteristic 

1 autocorrelation time 

AvgDif! » SumDif# / (Jfc - 1)* Average absolute value change in seq readings 
RmaDif] m SQR(SumDif2# / (J* - !))• rks change in sequential readings 


1 Next go through the data array a couple of times. 

• Find and characterize large events (crossings of average and 

• avg+ (std-dev) lines) . 

1 Note-, the following is not the float efficient method of performing this 
1 analysis, but is logically simple and easier to read (and debug) . A final 
1 version could be performed in half as many operations- 

• crossings of the average value line: 

• First, find the first positive crossing of average line: 

X - 1 
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Y - 0 

DO UNTIL Y - 1 

IF axray<X + 1 - 30000) > Av*j and array (X * 300001 < Ave! THEN T - l 
X • X +• 1 

LOOP 


’Re suiting X is the first positive crossing of the average line 

• Next, stare there sad find number of positive crossings and duration of 

• periods abovs average and periods between positive crossings 


aNuct 0 
avgCroea « 
ADur * 0 
aDur2& - Q 
aPd2fc - 0 


*?. +ve crossings counted yet 
Pointer to most recent *ve crossing of avg 
1 Accumulator for duration above avg 
* Accumulator for duration squared above avg 
1 Accumulator for square-length of crossing period 


FOR ptNum - X * 1 TO NArray 

IF array tpcNUn - 30000) » Avel AND array (ptSum - 1 . JOOOOI < Aval THEN 
* +ve exosaing 

aNum ■ aNum ♦ 1 * Count the crossing 

aPdlSc ■ aPd2A * (ptNum - avgCroea} * 2 'Accumulate square period 
avgCrosa - ptNum 'Update pointer to +ve crossing 
END IF 

IF array (ptNum - 30000 « Av«| and array (ptNum - l - 30 OOO) > Avel THEN 
' -ve c roe sing 

ADur - ADur ♦ ptNum - avgCrosa • Accumulate duration above bar 
aDur2* - aDur2t * (peNum - avgCrosa) * 2 * Aecum dor *2 above bar 

last Minus » ptNum 'store pointer to *v» crossing 
END IF ' 

NEXT ptNum 


» Calculate statistics for average-crossing events 

APdi * (laatMinus - X) / aNum • avg period of average -crossing events 

AKxisPd! * SQR(aPd2fc / aNum) * RMS • — 

ADuri - ADur / aNum 1 Avg duration of time above avg 

ARmaDur! - SQR(aDur2I / aKum) 4 sms — 


• Next, look for crossings of the average* ( standard deviation) line: 

* First , find the first positive crossing of avg+tstd-dev) line? 

X - 1 
Y - 0 

bar! *» Avel * SDevi 
DO UNTIL Y « 1 

IF array (X * 1 - 30000) > bar! AND array(X - 30000) < bar! THEN Y * 1 
X - X ♦ 1 X 

LOOP 

'Resulting X ic the first positive crossing of the avg*( 6 td-dev) line 

♦ Next, start there and find number of positive crossings and duration of 
1 periods above the bar and periods between positive crossings 

aNum - 0 * No +ve avg* (std-dev) crossings counted yet 

avg Croce - X • Pointer to most recent +ve crossing of avg* (std-dev) 

ADur • 0 1 Accumulator for duration above avg* (std-dev) 
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aDur2fc - 0 ’ Accvaailator for duration squared above *vg+ (std-dav) 

aPd2* - 0 * Accumulator for square* length of crossing period 

FOR ptNum * X ♦ 1 TO NArray 

IF array (ptNum - 30030) > bar l AND array (ptNum - 1 - 30000) < bar! THEM 
1 +ve crossing 

aNura « sMuro + 1 'Count the crossing 

aFd2& » *Pd2£ ♦ (ptNum - avgCross) * 2 'Accumulate square period 
ftvgCross - ptNoa 'Update pointer to +ve crossing 
END IF 

IF array (ptNum - 300QQ) < barl ASD array (ptNum * I - 30000) > bar! THEN 
1 -ve crossing 

ADur * ADur + ptNum - avgCzoea • Accumulate duration above bar 
aDurlA • aDur26 + (ptNum * avgCross) " 2 1 Aecum dur"2 above bar 

lastMinua ■ ptNum 'Store pointer to -ve crossing 
END IF 
NEXT ptNum 

* Calculate statistics for avg+ (atd-dev) crossing events 

SPdi • (last M inu s - X) / aNum ' Avg period of events above avg* (std-dev) 

SRmsFdl - SQR UPd2i / aNum) -rms ■ 

SDur ! - ADur / aNum *Avg duration of time above bar 

SRmsDur! « SQR(aDur2fc / aNum) f HMS — 1 

• Print results in output file: 

1 First, find a. 3 format name of input file Ci*«** atrip off path) 

filename $ « ** ' Preset to null in case path isn't present 

FOR i . 1 TO l*EN(f$ (filaNum) ) 

IF MIPS tf$ (fileNum) , LEN(f$(fUeiium)) ♦ 1 - i, l) * «\* THEH 
fileName$ - RIGHT? (£S CfileNum) , i - 1) 

EXIT FOR 
END IF 
NEXT i 

IF me»aine$ • — TEEN filename? - £$ (fileNum) 1 No path 
SELECT CASE LEN(£lleKam$): 

CASE 12: ' Sheesh, what we won’t do for A clean output format 1 

PRINT 42, fileName?; 

CASE Us 

PRINT 42, * ■; filename? ; 

CASE 10; 

PRINT 42, " ■; £ileNa*e$ ; 

CASE 9z 

PRINT 42 r * *; fiiaKameS; 

CASE 8; 

PRINT #2, * *; filename? ; 

CASE 7: 

PRINT 42, 0 •; fil«Name$ ; 

CASE Gt 

PRINT 42, ■ *s fUeSan»$. 

CASE 5 : 

PRINT # 2 , ■ m ; IUcNameS; 

CASE 4 1 

PRINT 42, * *? fileNaoeS- 

CASE 3; 

PRINT 42, " file«ame$ r 
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CASE 2: 

PRINT M2, - 

CASE Is 

PRINT #2. ■ *7 

END SELECT 

PRINT #2, USING "*###*.# #f ; Ave!; 

PRINT #2, • 

PRINT #2, USING "*«#44.ffti' ; SDev!j 
PRINT #2 , • •; 

PRINT #2, USING ■##*##.«#*; Antotf 
PRINT #2, • •? 

PRINT 12, USING *##!»•.*#*• ; AvgDif I j 
PRINT #2, ■ •; 

PRINT #2, USING "#*§#*.«•«'; taaDlfls 
PRINT #2, * *; 

PRINT *2, USING ’*##«#.*##"; APdl; 

PRINT #2, • •; 

PRINT #2 , USING *#####.###* i ARnsPdt ; 

PRINT #2. ■ •* 

PRINT #2, USING •Mfl#0.«ftt"i ADur 1 ; 

PRINT #2, • *s 

PRINT #2, USING ”#*«#! •*##*; ARnaDurl : 

PRINT #2, * *} 

PRINT #2, USING Wii Ui*; SPdl f 
PRINT #2. ■ *; 

PRINT #2, USING ■####!.«*•; SRasPdl ; 

PRINT #2, - *7 

PRINT #2, USING •#####-##*•; SDurlj 
PRINT #2, • "i 

PRINT M2, USING •#####.###*; SRmsDur! 

Restart: ‘Re-entry point Cor Kissing input file 

NEXT fileNum 
CLOSE #2 
END 

Handle! : ' Fatal errors that require program termination: 

PRINT ^ 

IF flag * 1 TB£» PRINT "Couldn't open list file *? FileListS; • for in 

IP flag - 2 THEN PRINT "No file names in FtleLiatS, 

IF flag • 1 THEN psun "Couldn't open •; oFiie$; • for output." 

PRINT 

END 

Handle2 : 1 Non- fatal errors that allow progr a m to continue : 

PRINT 

PRINT "Can't open file f$ ; * f or input." 

ON ERROR GOTO 0 
OOTO Restart 


defint A-Z 

SUB Comline {NumArgs, ArgB$ 0 , MaxArgs) STATIC 
CONST TRUE - -1, FALSE - 0 

NumArga » 0 j in - FALSE 
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> Get the command lisa using the OOMMAKD$ function. 
Cl$ - CQMKAKD$ 
l» m LEU (Cl$ J 

1 Go through the conwrTin d line a character at a time* 
FOR i - 1 TO L 

C$ - MXD$ (cl$, i, 1) 

‘Test for character being a blank or a tab. 

IF (C$ <> • • AMD C$ CHR$19) ) THEM 

* Neither blank nor tab. 

* Teat to see if you're already 

* inside an argument. 

IF ROT in THEN 

• You've found the start of a new argument. 

' Test for too many arguments. 

IF NumArga «• HaXArga THEN EXIT FOR 
NumArgs - NunAxgs + x 
in * TRUE 
END IF 

1 Add the character to the current argu me nt . 

Arg«$(Nun*rgs) « Args$<NmoArga) ♦ C$ 

ELSE 

* Found a blank or a tab. 

' Set *Not in an argument” flag to FALSE. 

In - FALSE 
END IF 
NEXT i 

END SUB 
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This program accepts a command line of the form: 
NAS96072 [DataFile$][-l FileList$]. If the analysis of a 
single data file is of interest then the DataFileS field above 
is filled with the path and name of the file containing the 
data. If a list of files are to be analyzed, then the -1 switch 
is used and the path of the name of an ASCII file containing 
a list of the files used as the FileList$ variable. In either 
event, single ASCII output files are created from the Data- 
FileS or ListFile$ name, with the file name extension 
replaced with “OUT”. This file contains one row that is a 
header describing the names of each of the signature quan- 
tities that was calculated, plus one row for each data file that 
was analyzed containing the values of the corresponding 
signature quantities. 

Neural network 80, FIG. 8, is then used to determine the 
flow rate and other typical variables associated with single 
and two phase flows in conduit 14, FIG. 1. The primary 
reason neural network 80 is used is that it has the ability to 
handle arbitrarily non-linear correlations involving abrupt 
transitions, strongly varying behaviors, or complicated inter- 
actions among input variables. In addition, neural network 
80 does not require any assumptions about the form or the 
correlation between a set of inputs and the desired set of 
outputs. Neural network 80 is programmed by training a set 
of weighted connections, rather than by explicitly defining a 
set of operations that must be performed. Once trained, 
neural network 80 is then implemented in a very fast, 
efficient software, firmware, or hardware environment ren- 
dering it ideally suited for real time applications. In addition, 
once properly trained, neural network 80 is extraordinarily 
fault tolerant, and/or fault indicating. Thus, if the data is 
noisy or occasionally erroneous or missing, a neural network 
can indicate this and still provide a reasonable estimate of 
the operating conditions. 

In practice, neural network 80 is trained by first inputting 
known flow indicator quantities for a broad range of known 
conditions and training the network to handle those particu- 
lar cases. Once network 80 is trained, it can be applied to a 
new set of flow indicator quantities providing a means of 
analyzing data in a real world applications. 

Neural network 80 is best viewed as being made up of a 
layer of neurons. The most common type of a neural 
network, called a “back propagation” network, contains 
three layers: input layer 82, hidden layer 84, and output layer 
86. The input layer receives the input information, in this 
case the flow indicator quantities discussed above. Each of 
the inputs is scaled from its full range to a range of (-1, 1,) 
which the network is very good at handling. No processing 
is done at the input layer, instead each neuron simply passes 
along its input value which is fed to the next layer through 
a set of weighted connections. 

Usually, successive layers are thoroughly interconnected 
so the weighted connections link each neuron in the input 
layer to each neuron in the hidden layer. The value that 
reaches a hidden neuron along a given connection is the 
product of the value injected into the connection by the input 
neuron and the value of the connection weight. In hidden 
layer 84, each neuron sums the values that reach it via the 
connections and then applies a simple, non-linear transfor- 
mation to that value. Gaussian (bell-shaped) and sigmoidal 
(S-shaped) functions are the most commonly used partly 
because they easily calculate. 

In a preferred embodiment, neural network 80 uses a 
sigmoidal function called a “logistic” function. This func- 
tion maps any input value in an infinite range to an output 
of (-1. 1). The resulting value is then fed forward along the 
weighted connections to output layer 86. Output layer 86 
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functions identically to the hidden layer, summing the 
weighted values that are presented to it, applying a logistic 
transformation, and outputting the result. The output value 
in the range (-1.1) is scaled back to the full range of interest, 
5 producing a final result of the network. 

In the subject invention, the eighteen flow indicator 
quantity inputs were used for each data set and one or two 
outputs corresponding to the flow rates or yes/no answers. 
The number of hidden neurons can vary considerably. 
10 depending upon the difficulty of the required calculation. 
More neurons permit greater complexity, but slow the pro- 
cessing speed. The “NeuroShell 2” neural network simula- 
tion program available from Ward Systems Group may be 
used to implement neural network 80, FIG. 8. 

15 Training back propagation network 80, requires the prepa- 

ration of a set of training data sets containing both the inputs 
and the desired outputs. The weighting values initially are 
selected at random. During training, training cases are 
chosen at random, presented to the network, and the error in 
20 the network’s output is “back propagated” through the 
network to determine minimum changes to the weighted 
values that will eliminate the error. By performing this 
process over and over, updating the weighted values as 
needed, the network’s calculational accuracy continually 
25 improves until, finally, it is not improved further. If the 
network’s accuracy is acceptable at that point, it can be used 
as desired to rapidly produce an output value for the set of 
input values. 

Neural network 80 may be subject to “over-training” 
30 unless care is taken to train it properly. Over-training is the 
ability of the network to recognize individual training sets, 
rather than the basic variation in the patterns between 
different training sets. Over-trained networks often perform 
poorly when used to estimate the state of new input sets. 
35 This problem can be avoided by using one sub-set of the 
available data sets as training data, and the remainder as test 
data. The network is trained using the training sets, but a new 
set of values for the connection weights is saved only when 
it performs better with the test data (which were not used to 
40 train the network). This method usually results in a network 
with a satisfactory ability to “generalize” beyond the train- 
ing data to the larger range of possible inputs. 

Each data file of raw data is filtered to a signature file 
using the computer program described above. These signa- 
45 ture files are merged with a master spreadsheet that con- 
tained the pressure, temperature, and flow rate information 
for each test. From this master spreadsheet, the appropriate 
information with respect to the phase content (liquid/two- 
phase), flow (flow/no-flow) and flow indicator quantities are 
so input to neural network 80 for processing. 

To train neural network 80, the water flow of a two phase 
water/air flow is determined using an electric digital flow 
meter such as a Great Plains (Wichita, Kans.) 0.3-3 gpm 
and/or 3 gpm-50 gpm capacity meter. Air flow is measured 
55 using either or both of two mass flow meters from Omega 
Engineering (Stamford. Conn.) with flow ranges of 0-50 
standard liters/minutes of nitrogen. 

For air- water flows, the following parameters are of 
primary interest: whether there is a flow in the pipe, whether 
60 there is air in the flow, if the flow is pure liquid, the flow rate: 
If the flow is two-phase, the water and air flow rates are also 
of primary interest Eighty-two test sets are selected at 
random as test sets, leaving 121 sets that may be used for 
training the network. Whether there is flow in the pipe and 
65 whether there is air in the flow are both simple yes/no 
answers. So, a value of 1 is set to the desired output for a true 
answer and 0 for a false answer. In actuality, any value 
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greater than 0.5 for a true answer or less than 0.5 for a false 12. The system of claim 9 in which a said flow indicator 
answer is effectively as good. quantity includes the root mean square value of the duration 

Once neural network 80 is trained in accordance with this of average crossing waves variable, 
invention, it will predict the existence of a flow and whether 13. The system of claim 4 in which a said flow indicator 


there was presence of air in the flow. 

Another neural network may then be trained, for pure 
liquid flows, to determine the water flow rate. 

Another neural network may be trained in accordance 
with this invention to determine the water and air flow rates 
of a two-phase flow. 

Thus, the subject invention reduces the sampling rate 
required by generating an amplitude envelope of the sensor 
signal thus eliminating the need for complicated hardware 
and software. And, the subject invention eliminates the need 
to use invasive sensors which often fail to work properly at 
low or high temperatures and in harsh or adverse conditions. 

Although specific features of this invention are shown in 
some drawings and not others, this is for convenience only 
as each feature may be combined with any or all of the other 
features in accordance with the invention. 

Other embodiments will occur to those skilled in the art 
and are within the following claims: 

What is claimed is: 

1. A flow analysis system comprising: 

a sensor in communication with a conduit including a 
flow therein for transmitting a signal which varies 
depending on the characteristics of the flow in the 
conduit; 

a filter, responsive to said signal, tuned to pass a narrow 
band of frequencies proximate the resonant frequency 
of the sensor; 

means for generating an amplitude envelope of the filtered 
signal; 

means, responsive to the amplitude envelope, for calcu- 
lating a plurality of flow indicator quantities represen- 
tative of the flow in the conduit; and 

means, responsive to the flow indicator quantities, for 
determining the flow rate of the flow in the conduit. 

2. The system of claim 1 in which said means for 
calculating includes means for sampling the amplitude enve- 
lope at a rate less than or equal to the resonant frequency of 
the sensor. 

3. The system of claim 1 in which said means for 
determining includes a neural network trained to calculate 
the flow rate in accordance with a plurality of flow indicator 
quantities representative of a known flow rate in the conduit. 

4. The system of claim 1 in which said means for 
calculating includes means for determining variations in the 
amplitude of the amplitude envelope. 

5. The system of claim 4 in which a said flow indicator 
quantity includes the average value of the amplitude. 

6. The system of claim 4 in which a said flow indicator 
quantity includes the standard deviation of the amplitude. 

7. The system of claim 4 in which a said flow indicator 
quantity includes the characteristic autocorrelation time. 

8. The system of claim 4 in which a said flow indicator 
quantity includes the average absolute value difference of 
the amplitude. 

9. The system of claim 4 in which a said flow indicator 
quantity includes the average crossing waves variable. 

10. The system of claim 9 in which a said flow indicator 
quantity includes the root mean square difference of the 
average crossing waves variable. 

11. The system of claim 9 in which a said flow indicator 
quantity includes the duration of the average crossing waves 
variable. 


5 quantity includes the period of standard deviation crossing 
waves. 

14. The system of claim 13 in which a said flow indicator 
quantity includes duration of crossing waves. 

15. The system of claim 4 in which a said flow indicator 
10 quantity includes the root mean square duration of standard 

deviation crossing waves. 

16. The system of claim 4 in which a said flow indicator 
quantity includes both the average value of the amplitude of 
the envelope and the standard deviation of the average value 

15 of the amplitude of the envelope. 

17. A flow analysis system comprising: 

an acoustic sensor in communication with a conduit 
including a flow therein for transmitting a signal which 
varies depending on the characteristics of the flow in 
20 the conduit; 

means, responsive to the amplitude envelope, for calcu- 
lating a plurality of flow indicator quantities represen- 
tative of the flow in the conduit; and 

means, responsive to the flow indicator quantities, for 
25 determining the flow rate of the flow in the conduit 

18. The system of claim 17 further including an amplifier 
for amplifying said signal. 

19. The system of claim 17 further including a filter, tuned 
to pass a narrow band of frequencies proximate the resonant 

30 frequency of the sensor. 

20. The system of claim 17 further including means for 
generating an amplitude envelope of the filtered signal. 

21. The system of claim 20 further including means for 
sampling said amplitude envelope at a rate less than or equal 

35 to the frequency of the signal transmitted by the sensor. 

22. A flow analysis system comprising: 

an acoustic sensor in communication with a conduit 
including a flow therein for transmitting a signal which 
varies depending on the characteristics of the flow in 
40 the conduit; 

means for amplifying the signal; 

a filter, responsive to the amplified signal and tuned to 
pass a narrow band of frequencies the proximate reso- 
45 nant frequency of the acoustic sensor; 

means for generating an amplitude envelope of the filtered 
signal; 

means, responsive to the amplitude envelope, for calcu- 
lating a plurality of flow indicator quantities represen- 
50 tative of the flow in the conduit; and 

means, responsive to the flow indicator quantities, for 
determining the flow rate of the flow in the conduit. 

23. The system of claim 22 in which said means for 
calculating includes means for sampling the amplitude enve- 

55 lope at a rate less than or equal to the resonant frequency of 
the sensor. 

24. The system of claim 22 in which said means for 
determining includes a neural network trained to calculate 
the flow rate in accordance with a plurality of flow indicator 

60 quantities representative of a known flow rate in the conduit. 

25. The system of claim 22 in which said means for 
calculating includes means for determining variations in the 
amplitude of the amplitude envelope. 

26. A flow analysis method comprising: 

65 coupling a sensor to a conduit with a flow therein and 
transmitting a signal which varies depending on the 
characteristics of the flow in the conduit; 
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using a filter, responsive to the signal and tuned to pass a 
narrow band of frequencies proximate the resonant 
frequency of the sensor; 

generating an amplitude envelope of the filtered signal; 

calculating, in response to the amplitude envelope, a 
plurality of flow indicator quantities representative of 
the flow in the conduit; and 

in response to the flow indicator quantities, determining 
the flow rate of the flow in the conduit. 

27. The method of claim 26 in which calculating includes 
sampling the amplitude envelope at a rate less than or equal 
to the resonant frequency of the sensor. 

28. The method of claim 26 in which determining 
includes using a neural network trained to calculate the flow 
rate in accordance with a plurality of flow indicator quan- 
tities representative of a known flow rate in the conduit. 

29. The method of claim 26 in which calculating includes 
determining variations in the amplitude of the amplitude 
envelope. 

30. The method of claim 29 in which a said flow indicator 
quantity includes the average value of the amplitude. 

31. The method of claim 29 in which a said flow indicator 
quantity includes the standard deviation of the amplitude. 

32. The method of claim 29 in which a said flow indicator 
quantity includes the characteristic auto-correlation time. 

33. The method of claim 29 in which a said flow indicator 
quantity includes the average absolute value difference of 
the amplitude. 

34. The method of claim 29 in which a said flow indicator 
quantity includes the average crossing waves variable. 

35. The method of claim 29 in which a said flow indicator 
quantity includes the root mean square difference of the 
average crossing waves variable. 

36. The method of claim 29 in which a said flow indicator 
quantity includes the duration of the average crossing waves 
variable. 

37. The method of claim 29 in which a said flow indicator 
quantity includes the root mean square value of the duration 
of average crossing waves variable. 

38. The method of claim 29 in which a said flow indicator 
quantity includes the period of standard deviation crossings. 

39. The method of claim 29 in which a said flow indicator 
quantity includes the duration of crossing waves. 

40. The method of claim 29 in which a said flow indicator 
quantity includes the root mean square duration of standard 
deviation crossing waves. 

41. The method of claim 29 in which a said flow indicator 
quantity includes both the average value of the amplitude of 
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the envelope and the standard deviation of the average value 
of the amplitude of the envelope. 

42. A flow analysis method comprising; 

coupling an acoustic sensor to a conduit with a flow 
therein and transmitting a signal which varies depend- 
ing on the characteristics of the flow in the conduit; 

using a filter, responsive to the signal and tuned to pass a 
narrow band of frequencies proximate the resonant 
frequency of the acoustic sensor; 

generating an amplitude envelope of the filtered signal; 

calculating a plurality of flow indicator quantities repre- 
sentative of the flow in the conduit from the amplitude 
envelope; and 

in response to the flow indicator quantities, determining 
the flow rate of the flow in the conduit. 

43. The method of claim 42 in which calculating includes 
sampling the amplitude envelope at a rate less than or equal 
to the resonant frequency of the sensor. 

44. The method of claim 42 in which determining 
includes using a neural network trained to calculate the flow 
rate in accordance with a plurality of flow indicator quan- 
tities representative of a known flow rate in the conduit. 

45. The method of claim 42 in which calculating includes 
determining variations in the amplitude of the amplitude 
envelope. 

46. A flow analysis method comprising: 

coupling an acoustic sensor to a conduit with a flow 
therein and transmitting a signal which varies depend- 
ing on the characteristics of the flow in the conduit; 

using the amplitude envelope to calculate a plurality of 
flow indicator quantities representative of the flow in 
the conduit; and 

determining the flow rate of the flow in the conduit from 
the flow indicator quantities. 

47. The method of claim 46 further including amplifying 
the acoustic sensor signal. 

48. The method of claim 46 further including using a 
filter, tuned to pass frequencies proximate the resonant 
frequency of the acoustic sensor. 

49. The method of claim 46 further including generating 
an amplitude envelope of the acoustic sensor signal. 

50. The method of claim 49 further including sampling 
the amplitude envelope at a rate less than or equal to the 
resonant frequency of the acoustic sensor. 

* * * * * 
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